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We study the effect of superconducting fluctuations on the longitudinal and the transverse (Hall) 
conductivity in homogeneously disordered films. Our calculation is based on the Usadel equation in 
the real-time formulation. We adjust this approach to derive analytic expressions for the fluctuation 
corrections in the entire metallic part of the temperature-magnetic field phase diagram, including 
the effects of both classical and quantum fluctuations. This method allows to obtain fluctuation cor- 
rections in a compact and effective way, establishing a direct connection between phenomenological 
and microscopic calculations. 



PACS numbers: 74.25.Fy, 74.76.-w, 74.40.-n 

I. INTRODUCTION 

Theoretical studies of fluctuation conductivity in su- 
perconductors found their origin in the discovery of para- 
conductivity by Aslamazov and Larkin (AL) in 196^. 
These authors analyzed the conductivity of supercon- 
ductors in the metallic phase above the transition tem- 
perature Tc in the framework of diagrammatic linear re- 
sponse theory. Paraconductivity can be understood as 
the direct contribution of fluctuating Cooper pairs to the 
electric current. Indeed, the formation of Cooper pairs 
opens a new channel for charge transport in the metallic 
phase. Above the transition temperature, these Cooper 
pairs do not form a condensate yet and their contribu- 
tion to conductivity is positive but still bounded due to 
their finite lifetime. Other effects of superconducting 
fiuctuations are Andreev scattering of electrons off the 
fluctuating order parameter de scribed by the so-called 
Maki-Thompson (MT) ternP^, and the suppression of 
the quasiparticle density of states (DOS) near the Fermi- 
level. 

These classical results were obtained for temperatures 
close to Tc and later extended for larger temperatures 
and for weak magnetic fields. More recently, the vicinity 
of the magnetic field-tuned quantum phase transition in 
disordered superconducting films was studied in a paper 
by Galitski and Larkirp'. These authors have shown that 
close to the quantum transition, contrary to the previ- 
ously studied regime of weak magnetic fields, different 
processes are of equal importance. This has the remark- 
able consequence that the sign of the total correction to 
conductivity becomes negative for sufficiently low tem- 
peratures near the quantum critical point, resulting in a 
non-monotonic magnetoresistance in this regime. 

In spite of the substantial amount of existing theoreti- 
cal work on superconducting fluctuations, summarized in 
the book by Larkin and Varlamov^, the subject continues 
to be an active field of research. This activity is stimu- 
lated by recent accurate experimental studies of different 



superconducting systemJ^H^i', that call for refined theo- 
retical studies. For example, when fitting experimental 
data on disordered superconducting films by theoretical 
results, one commonly uses several fitting parameters, 
including the critical temperature Tc, the upper critical 
field Be and the dephasing time r^. In doing so, it would 
be useful to work with theoretical results which are valid 
in the entire {B, T) phase diagram, instead of address- 
ing different asymptotic regions separately. This is the 
motivation for the detailed calculations presented in this 
paper. 

In deriving the results for the fluctuation conductivity, 
we deviate from the traditional route that employs the di- 
agrammatic linear response theory in the imaginary time 
techniqu^i^l as described in detail, for example, in Ref. [S] 
Instead, we develop a formalism based on the Keldysh 
(real-time) representation of the Usadel equation. In this 
approach, disorder averaging is performed at the earliest 
stages, thereby avoiding the use of the impurity-diagram 
technique. As an additional advantage, no analytic con- 
tinuation is required. The Usadel equatioiP^ is an indis- 
pensable tool in the theo ry of mesoscopic superconduc- 
tors and hybrid structure^^E^. This equation describes 
low-energy (diffusive) physics on spatial q^^ and tem- 
poral uj~^ scales, satisfying {ql, ujt) ^ 1, where r is 
the impurity scattering time and I the mean free path. 
The first calculation of superconducting fluctuation cor- 
rections in this framework was performed by Volkov et 
al^, who calculated fluctuation conductivity in hybrid 
superconducting-normal structures in the vicinity of T^. 
in the absence of a magnetic field. 

In this paper, we use the Usadel equation to calculate 
longitudinal and transverse (Hall) conductivity in disor- 
dered superconducting films at arbitrary temperatures 
and magnetic fields. Our approach parallels to some ex- 
tent the non-linear cr-model formalism for disordered su- 
perconductors introduced by Feigelman et al.'^, and the 
subsequent work by Kamenev and LevchenkcP^. The lat- 
ter work includes a calculation of fluctuation conductivity 
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close to Tc- The intimate relation between the cr- model 
formalism and the Usadel equation approach is based on 
the fact that the Usadel equation is the saddle point equa- 
tion of the nonlinear cr-model. For the sake of simplicity, 
we decided not to use the more technical apparatus of 
the nonlinear a-model, but formulate the derivation in 
terms of the Usadel equation. This route leads us to a 
description in terms of a coupled set of kinetic equations 
for quasiparticles moving on the background of supercon- 
ducting fluctuations. This method appears to be a very 
convenient tool for studying fluctuation transport. 

The classification of the fluctuation corrections ob- 
tained in the discussed method appears to be very dif- 
ferent from the conventional classiflcation based on the 
diagrams in the Matsubara technique. Therefore com- 
parison with the results obtained by the diagrammatic 
technique can be performed only on the level of the flnal 
results. Let us mention here the comparison to recent 
works. It can be seerP^that the zero magnetic field limit 
of the general formulas derived in this manuscript (Eqs. 
(771, (78), and ([79| below) can be presented in a form 



that exactly coincides with the corresponding diagram- 
matic results of Lopatin et al. in Refs. I20|21l On the 
other hand, Glatz et al. more recently presented a di- 
agrammatic analysis of the longitu dinal fluctuation con- 
ductivity in the entire phase diagrarrp2123| However, their 
results are inconsistent with previous diagrammatic cal- 
culations as well as with ours (we comment on this work 
below at the end of Sec. IV B). For the Hall effect, our 
results agree with those of a worl?!^ in which an inde- 
pendent calculation has been performed. These results 
were successfully applied for the description of a recent 
measurement in amorphous Tantalum Nitride fllms.l^ 

This paper is organized as follows. In Sec.[lT]we present 
the basic formalism. We show how the Usadel equa- 
tion, initiall y f ormulated for a given order parameter 
configuration^, can be applied to the calculation of fiuc- 
tuation conductivity. As a next step, in Sec. |III| we find 
a solution of the Usadel equation which allows to de- 
termine the order parameter correlation function in the 
Gaussian approximation. Both ingredients are required 
for the calculation of the electric current presented in 
Sec. IV A Next, we derive expressions for the longitu- 
dinal conductivity that are valid in the entire metallic 
phase outside the regime of strong fluctuations. Eval- 
uation of the obtained expressions still requires a sum- 
mation over the Landau levels as well as an integration 
over slow (bosonic) frequencies, which can be performed 
analytically only in certain limiting cases. Several such 
limiting cases are analyzed in detail in Sec. |IVB[ includ- 
ing the region close to Tc and the vicinity of the quantum 
critical point. By means of a numerical evaluation, we 
locate the line of the sign change for magnetoresistance 
da/dB and the fine da/dT = 0. We also discuss the exis- 
tence of a crossing point of the magnetoresistance curves. 
In Sec. |IV C| we calculate Hall conductivity, generalizing 
previous calculation^SMUl to the case of arbitrary tem- 
peratures and magnetic fields above the transition. 



II. BASIC EQUATIONS 

In this section we present the equations that form the 
basis for our calculation of the fluctuation conductiv- 
ity. After stating the microscopic model, we introduce 
the Usadel equation that allows to find the quasiclassical 
Green's function in the dirty limit, i.e., if the condition 
TcT ^ 1 is fulfilled. Calculation of the conductivity re- 
quires knowledge of both the quasiclassical Green's func- 
tion in the presence of the fluctuating order parameter 
fleld and the correlation function of the order parame- 
ter field. In the fluctuation regime, which we study in 
this paper, the order parameter correlation function is 
governed by the Ginzburg-Landau (GL) action. Fortu- 
nately, the GL action can be found from the quasiclassical 
Green's function itself, i.e., from the solution of the Us- 
adel equation. This procedure will also be described in 
this section. 

We start with the Keldysh action for electrons with 
short-range BCS-type interaction. After decoupling the 
interaction with the help of a Hubbard-Stratonovich 
transformation, the resulting action is split into two parts 
5[*,A] = S^ii*, A] + S'2[A], where 

5i[*,A]= (1) 
j dx *^(x) [it:idt - H{x) + ^ + A {x)] *(a;), 
2v 



S2[K] = j dxtr [A+aiA] . 



(2) 



Here, v is the density of states per one spin projection 
at the Fermi level and fi is the chemical potential. The 
dimensionless coupling constant in the Cooper channel 
A is positive for an attractive interaction. Hereafter, we 
use the hat symbol as in to denote 2x2 matrices 
in Keldysh {K, retarded/advanced) or Gor'kov-Nambu 
{N, particle/hole) spaces. By ct^ and we denote the 
Pauli matrices in K and N space, correspondingly. The 
check symbol as in H denotes 4x4 matrices in the direct 
product space K ® N . The trace operation tr in Eq. ^ 
comprises both K and N spaces. The short notation 
X — (r, t) is used, and the time integration covers the 
interval (— c)o,oo). The single-particle Hamiltonian H is 
defined as 

= --^ (V - ieA{v)hf + [/ (r) + e(^(r), (3) 

with a static disorder potential U , scalar ip and vector 
potentials A, and electron mass m and charge e. In the 
action, ^ is a four component vector of Grassmann fields 
with the following structure: 



^1 
V'2 



4 A 



K 



K 



N 



V'l = (Xit,-X4) 



N- 



(4) 



(5) 
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All terms in the electronic action Si are diagonal in K- 
space except the order parameter field A = Aoao + AiCTi, 
where Aq and Ai will be referred as classical (cl) and 
quantum (q) components of the order parameter. These 
components are non-diagonal in N space: A^ = Ajf+ — 
A*f_, where t± — {t^ ± iTy). We arrange the classical 
and quantum order parameter fields into the vector A = 

(Ae/,A,r. 

The electronic Green's function for the system reads: 

This expression can be cast in the form 

Gix,x') ^ J DA Ga{x,x') e'^'^^I^l, (7) 
where the Ginzburg-Landau action is determined by 

5GL[A] = -zlnyi?*e'^[*'^l, (8) 

while 

GA{x,x) = -^ • (9) 

This Green's function depends on the specific configura- 
tion of the order parameter field A. 

Physical quantities can be obtained in terms of the 
disorder-averaged {^G{x,x')')^.^, which can be found as 

(G(x,x')),, = /i^A(G^(.,.')L. e'<^-[^l)-. 

(10) 

Here, we average the electronic Green's function sepa- 
rately from the bosonic action. This is a valid approxi- 
mation for films with dimensionless conductance g ^ 1; 
taking into account cross-correlations between the two 
terms would give corrections to the Drude conductivity 
of the order of 1/5^, while we are only interested in cor- 
rections of the order of 1/g. 

The electric current is related to the Keldysh compo- 
nent of {G{x,x'))^.^: 



e 
2m 



(Vr-V.)(G^(a:,x')),__,-^A, (11) 



where n stays for the density of electrons. 

In the following, it will be convenient to use the quasi- 
classical approximatiorPsHsU. The quasiclassical Green's 
function can be introduced as follows. First, one per- 
forms the Wigner transform of the disorder-averaged 
Green's function as 



energy variable ^ = ^ — fi which describes the distance 
from the Fermi surface: 



i 

I" .'-00 



(Ga (n {pF + ^/vf) , r, ^1,^2))^,, .(13) 



In this equation, vf denotes the Fermi velocity. 

In the diffusive limit higher angular harmonics are 
suppressed and a formulation in terms of the angular- 
averaged Green's function is possible: 



g{r,ti,t2)^ J dn g^{r,ti,t2) 



(14) 



The function g satisfies the nonlinear Usadel 
equatioiP^E^ 

DV U ■ Vg) - {f^du g} + i[A- eif, g] = 0, (15) 



where the symbol • is used to denote a convolution in 
time, i.e., integration in the intermediate time variable. 
The spatial derivative has the following structure: Vg = 
Vg — ie [T3A, g] . An important constraint imposed on the 
quasiclassical Green's function is that it has to satisfy the 
normalization condition 



{g-g){t,t') = iS{t-t'). 



(16) 



In what follows we are interested in Gaussian fluctua- 
tions. This means, that the film is considered to be not 
too close to the superconducting transition. The width of 
the non-Gaussian region is determined by the Ginzburg 
number Gi; in the case of disordered films Gi ^ g~^. 
The precise criterion for the range of validity of this ap- 
proximation depends on the quantity in question. Con- 
cerning transport phenomena, the non-Gaussian region 
is wider than for thermodynamics and has been esti- 
mated to be of the order of V Gi for the thermal phase 
transitiorPsl, i.e. it covers the temperature regime for 
which \T - Tel < y/GiTc. To the best of our knowledge, 
there is no such calculation for the quantum transition 
(a study of the effect of fluctuations on the critical mag- 
netic field exists^^. In this paper, we assume that we are 
always outside the region of non-Gaussian fluctuations. 

Let us now turn to the discussion of the Ginzburg- 
Landau action. As long as we are interested in Gaussian 
fluctuations, we need to know S'(3l[A] only up to the 
second order in A. Noting the relation 



5{Sgl[A] 



dis 



5A*ixi) 



itY [CT,f_(GA(a;i,a;i))^J 



2z/ 



(o-iA(xi))j 



(17) 



we can obtain 



(JO) \ / dis J 



(12) 

where r = ^(ri -I- r2), p = (ri — r2). Next, the quasiclas- 
sical Green's function is obtained by integration over the 



-^o-i4i,a;2 -I- n(a;i,X2) 



(18) 

A(X2), 
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where 



Ilij{xi,X2) = i 



disi 



(19) 



A=0 



Importantly, the appearing Green's function at coincid- 
ing times and space points is related to the quasiclassical 
Green's function, and we can write 



II*^(xi,X2) = TTiy 



6A,ix2) 



(20) 



A=0 



This result shows that knowledge of the quasiclassical 
Green's function, i.e., the solution of the Usadel equa- 
tion, also allows finding the GL action. This observation 
considerably simplifies the scheme of calculation of the 
Gaussian corrections. With the help of the GL action, 
in turn, one can obtain the order parameter correlation 
function, which is needed for the calculation of the cur- 
rent. 

The charge density and electric current are expressed 
in terms of the angular-averaged Green's function g in 
the following wajff^ 

p (r, t) = -eiy (2e(l> + |tr {a^g (r, t, i))) (21) 



and 



CKVD 



tr(f3<Tij (r,i,t)) , 



(22) 



with ] ~ g • Vg. The angular brackets in this equa- 
tions symbolize averaging with the action Sgl- Relation 
(22) follows from Eq. (11) in the diffusion approxima- 



tion. Aiming for the needed accuracy (the leading order 
approximation in it is sufficient to determine j up 

to the second order in the fluctuating field A before the 
expansion in the electric field is performed. 



III. SOLUTION OF THE USADEL EQUATION 
AND THE ORDER PARAMETER 
CORRELATION FUNCTION 



For practical calculations, one needs to resolve the nor- 



malization condition (16) for the quasiclassical Green's 
function explicitly. In the framework of a mean-field 
treatment, one works with the classical order parameter 
field Ac/ only. In this case (Ag = 0) the Green's function 
can be parameterized as 



9 = 



g^ 



with 5^ = g^-^-h-h-g^ and {g^ ■ g^%^t' = (g^ ■ g'^)t,t' = 
I6t-t'- However, in the presence of the quantum order 
parameter fluctuation (i.e., for finite A^) this structure 
is broken and a more general parametrization needs to 



be considered. In that case, one can generalize (23) to 



take into account fluctuations up to the second order in 
A: 

g^-h-g^ g^ -h-h- g^ -h- g^ -h- g^ 
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g^+f-h 



In particula r, the lower left corner of this matrix is not 
equal to zercpSMl. This parametrization has the following 
property: 



{9^y^{g^)=i5,^,+0{A') (25) 



The matrix 



he 
hh 



(26) 



is called generalized distribution function'^^. Matrices 
gZ,w (diagonal and appear only in the second order 
in A. This holds provided the distribution function h 
satisfies the following normal metal diffusion equation: 



dt + ieipTs, h — 



(27) 



For the purpose of our calculation, we may assume 
g^ = (f^ = 0. In the case of g^ the reason is the fol- 
lowing. For the calculation of the current, the Green's 
function needs to be inserted into the corresponding ex- 
pression ( 22 ) and subsequently averaged over order pa- 
rameter configurations. There can be two kinds of con- 
tributions to the current originating from g^ . First, if 
it is not combined with any other term arising due to 
fiuctuations, it should be averaged by itself. Since the 
lower-left corner of the averaged Green's function must 
equal zero in the Keldysh formalism {g^) = 0, contribu- 
tions of this first type vanish automatically. The second 
kind of contribution appears when combining g^ with 



other terms arising due to fluctuations in formula (22). 
Since g^ itself is already quadratic in A, this procedure 
generates contributions to the current which are at least 
of the fourth order in A. These terms are beyond the ac- 
curacy of our calculation. The same argument applies to 
contributions originating from g^ , only in this case the 
average (g^) does not vanish identically, but is O(E^), 



as discussed in Appendix |Aj Therefore, there is no need 
to keep g^ when studying linear response in the elec- 
tric field. To conclude, for the purpose of our calculation 
we may work with the simple parametrization given in 
Eq.[23| 

In what follows, we consider static and homogeneous 
electric E and magnetic B fields and find it conve- 
nient to work in a gauge with time-independent elec- 
tromagnetic potentials: E = — V0 and B = curlA with 
(23) (f) — — Er, A — (0,-82:, 0). Under these conditions and 
in the absence of superconducting fluctuations, the re- 
tarded and advanced sectors of the quasiclassical Green's 
function are diagonal in N-space and take a particularly 
simple form 



g'^ih,t2)^-g^{h,t2)^f35t,. 



(28) 
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For a closed system, i.e. in the absence of a connection 
to an external bath, the distribution function h is not 



yet fixed. Indeed, equation ( 27 1 has infinitely many solu- 



tions. In the presence of interactions, it is convenient to 
work with the distribution function corresponding to the 
state of local thermal equilibrium with spatially varying 
chemical potential: 



h = 



he 








K^h^'H{eTe4>{x)) (29) 



where 



■H(e) = tanh 



2T' 



(30) 



This particular choice is especially convenient for linear 
response studies, because deviations of {g^) from zero 
which arise due to interactions are pushed to the second 
order in the electric field. This considerably simplifies 
perturbation theory. Note that temperature is still arbi- 
trary and is determined by the heat balance with a sub- 
strate or with contacts. Meanwhile, by neglecting (5^) 
we dismiss the heating effect of the electric field. 

In the presence of superconducting fiuctuations, the 
quasiclassical Green's function acquires off-diagonal com- 
ponents in iV-space. For the analysis of the Gaussian 
fluctuation regime, the deviations from the simple form 
given in Eq. [28] are small and may be treated as a per- 
turbation. With this in mind, we resolve the remaining 
constraints (25) as: 



1-5/-/* 



5^ 



-i + i./-/* 
-r 1 



hi* ■ I 



(31) 
(32) 



From the solution of the Usadel equation it will follow 
that /, / etc. are C(A). The functions / and /* as 
well as / and /* are considered as independent: they 
become complex conjugates of each other only when 
is neglected. 



We introduce parametrization (24) into Eq. (15) and 



neglect terms of the third order in A. As a result, we 
find for / the equation C~^/ — V, where the operator 
C""'^ is given by 



(33) 



and the gauge invariant derivative is: V/ — 
(V — 2ieA) f. As one may notice, this equation describes 
the response of the field / to the order parameter A, 
which enters this equation in the following combination: 



Vt,,t2 (r) = 2i [Aci (r, h) 5t,-u + K (r, h - ^2) A, (r, iz)] • 

(34) 

Similar equations arise for /*, /, and / with appropri- 
ately modified differential operators and functions V* , V 



and V . Taking into account the explicit form of he,h 
one may conclude that /ti,t2 = — /t2,ti (ttis same prop- 
erty holds for /*). Note that a static electric potential 
does not enter the equation for /. This is one of the 
advantages of the gauge in which the electric field is ex- 
pressed through the scalar potential. 

The equation for / can easily be solved after a Fourier 
transformation to the frequency domain according to 

/ (^1,^2) = / / (ei, £2) e-^('i*^-^^*^) (dei) (A2) . (35) 



Here we used notation (de) = de/2'n. To proceed further, 
we pass to the Landau level (LL) basis with eigenfunc- 
tions V'np (r) of the kinetic energy operator 



-D{V - 2ieA) V«p (r) = e„Vnp (r) 



(36) 



This equation describes a " particle" with a mass equal to 
1 /2D. We choose to work in the Landau gauge, for which 
the eigenfunctions ipnp are numbered by the momentum 
p and LL number n: 



(37) 



with magnetic length Ib = ^/ 
with charge 2|e|) and 



/2\e\B (for a "particle" 



Xn {x) 



1 



(x/Ib) ■ 



(38) 



Note that a description based on the Usadel equation is 
valid as long as we consider the regime of classically weak 
magnetic fields, for which ujc = satisfies WcT ^ 1. 
This means that the quantization of the orbital motion of 
the quasiparticles can be neglected. In contrast, the LL 
quantization of the collective modes and Cooperons e„ — 
f2^ (i 4- with ric — A\e\DB may still be important in 
the region of magnetic fields and temperatures we are 
interested in. 

The solution for / is conveniently written in terms of 
the Cooperon propagator, which is diagonal in the chosen 
basis: {n,p\C \m,p) = (5m„C„(ei + £2) with 



(39) 



Here, we introduced the dephasing time r^. The role of 
is to provide the long-time decay of the Cooperon, which 
is necessary to render corrections due to single-particle 
interference processes finite. These processes include 
weak localization and the anomalous Maki- Thompson 
correction (an analog of weak antilocalization) that di- 
verge in the absence of a magnetic field for r^^ =0. De- 
phasing can be provided by magnetic impurities or inelas- 
tic processes, i.e. electron-electron or electron-phonon 
collisions. For low temperatures, electron-electron colli- 
sions dominate. Outside the region of strong fluctuations 
(i.e., in the Gaussian regime), one can consider the de- 
phasing rate as energy independent and equal to the sum 
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of rates due to the CoulomtPSI and Cooper channel^^IISll 
In our study, we do not specify the dominant dephasing 
mechanism relevant for and consider it as an indepen- 
dent parameter. 

The solution of the equation C^^ f — V ioi f reads: 



and 



fnp (ei, £2) = C„ (2e) J Fe^.e, (r^Cp (r') dr' , (40) 
where 

.es (r) = 2i [Aci (r , w) + /le (r , e + w/2) A, (r , w)] 

(41) 

with shorthand notation e — (£1+62) /2 and w — ei — £2- 
Analogous equations hold for /*, / and /*. 

Having found approximate solutions for and , we 
turn to the GL action SVjl • As follows from Eq. [18] in 
combination with Eq. 20, it is sufficient to know ^^^^ 
at the first order in A to determine Sgl in the Gaussian 
approximation. We write the GL action in the form: 

Sgl[A] - J tr(2i.A+(~a.,r)/:-i(u;, r, r')A{uj,r') 



with 



£-1 



C^,' 

'-21 '-22 



(42) 



(43) 



Arguments (w, r, r') of £~"'^ are omitted in what follows. 

A straightforward calculation according to Eq. ( 18 ) 
gives: 



^21 i^np (r) Tplp (r') 



-w/2+e0(r) 



de 



1 



2e + i(e„+T, 1) A 



(44) 

The rest of the elements of £ ^ are related to this one 
according to = {'^21) ^^'^ 

C^i =B{u:-e4> (r) - e0 (r')) [C^l - C^^] , (45) 

where the bosonic distribution function is defined as 

S(w) = coth(a;/2r). (46) 

One can see, that the components of C^^ are not inde- 
pendent. Just as the components of £, they are related 
by the generalized fluctuation-dissipation theorem, see 
Eq. (45), valid in the quasi-cquilibrium state. Thus, only 
C^^ needs to be calculated explicitly. The evaluation of 
the e integral in Eq. (44) is straightforward and yields: 



= E V'«p(r)V:p(r')^n(w - 2e0(r)), 



where 



f„(u;)=ln^+V('^ ) -^^(n,c.) 



(47) 



(48) 



(49) 



We have introduced the BCS transition temperature 
Tc — ^^*^exp(— i), where ujd is the Debye frequency 
and 7 « 1.78. The symbol ip stands for the Digamma 
functioiP^. In deriving asymptotic expressions, we will 
use the following properties of this function: 2p'{l/2) = 
tt'^/2 and ip{x) sa Inx for x ^ 1. 

The line of the superconducting transition on the mean 
field level is determined by the equation £n=oi^ = 0) = 0. 
In the absence of dephasing — 00 this gives for the 
upper critical field 



BciT = 0) 



2^' 



(50) 



Let us discuss the effect of dephasing on the transition 
line. Since the fiuctuation propagator depends on the de- 
phasing rate, the transition temperature is shifted due to 
T0. Furthermore, since depends on the magnetic field 
as well as on the temperature, the presence of in the 
fiuctuation propagator changes the shape of the transi- 
tion line as a whole. Dephasing also affects the magne- 
toconductivity. This effect has been taken into account 
in the analysis of the experimental data on magnetocon- 
ductivity of thin superconducting InO films^^. 

As can be seen from the right-hand side of Eq. ( 47 ) , 



C21 is not translation invariant. However, by splitting 
off a gauge-dependent factor it can be rewritten in the 
following form: 



(51) 



/:2-i'(i,r,r') = e-'^'(*'-^-"-')£2"i'(^,r-r'), 
where Sg is defined as 

Sg{t,r,r') = e(0(r) + </)(r'))t 

-e(A(r)-KA(r'))(r-r') (52) 

and £-21 is translational and gauge invariant. We never- 
theless prefer to work with the operator in its origi- 
nal form. 

In order to find the order parameter correlation func- 
tions, one has to invert the operator C^^ given by 
Eq. (43) with the following result: 



C 



£^ 



where 

= (£2-/)-^ , £^ = (£r2^)-^ , = -£^£2-2^£-^. 

(54) 

The order parameter correlation functions are given by: 



K 



(Aez {io)A: (-^)) = ^/: 

(A,, (..)A; (-c.)) = ^£«, 

(A, (^)A:, i-u^)) = ^C^, 
(A, (^)A; (-..)) =0. 



(55) 
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In equilibrium, -C^^^q i'^) = L^^^^ (w) is diagonal in the 
LL basis, and reads as follows 



(56) 



For the Keldysh propagator this gives, according to 
Eq. (ImI): 



(u) = B (a.) (L« (c) - (..)) ^ (c.) . (57) 

While we have already neglected the heating induced 
by the electric field, we still keep other nonlinear effects. 
For example, one may consider the decay of fluctuat- 
ing Cooper pairs due to the acceleration of the paired 
electrons caused by the electric field. It was considered 
before on the basis of the phenomenological theorjffllEll 
(with only AL process included). At T ^ Tc this ef- 
fect becomes essential at electric fields of the order of 
E ~ Tc/e^GL that can be rather small due to the diver- 
gence of the coherence length ^g/, at the transition. 

In the following calculations all nonlinear effects will 
be neglected. In the linear response regime, we need to 
find propagators at first order in the electric field. Con- 
cerning the dependence of C on spatial arguments, we 
will consider it as an operator in the basis of the LLs, 
the same is assumed regarding the position operator r. 
Hence, in the equations below these two operators do 
not commute. We linearize £-21^ looking for the first or- 
der correction to its equilibrium value. In the equations 
below we do not indicate the frequency dependence of 
propagators, having in mind that all functions have the 
argument uj. Taking into account first-order corrections 
in the electric field we write 



= {l + 2eErd^)£. 

For £^ this gives: 

SL^ = ~2eEL^vdu:EL^ 



(58) 



(59) 
(60) 



and can be found by hermitian conjugation. Let us 
turn to . In order to find it, we need C22 given by 
Eq. p5l): 



C 



22 



B{C2^ 



^12 ) 



eEa^,S{(£-r),r}, (61) 



where curly brackets denote an anticommutator. Plug- 
ging this expression into Eq. ( 54 1 , we obtain 



-eEd,,BL^{{£ -£*),r} L^. 



(62) 
(63) 



Now, the order parameter correlation functions given by 
Eqs. (551 are fully specified, and we can proceed to the 



calculation of the electric current. 

To summarize, we have collected the basic elements 
of the formalism used for the calculation of the fluctu- 
ation conductivity in this paper. Once the quasiclassi- 
cal Green's function is found as a solution of the Usadel 



equation ( 15 ) , the current can be obtained from Eq. ( 22 1 



Since the quasiclassical Green's function is a functional 



of the order parameter configuration, formula (22 1 for 



the current includes an average with respect to the GL 
action. This action, in turn, can be found from the 
quasiclassical Green's function via Eqs. (18) and (20) 
and thus a closed scheme is established. ^As we have 
already argued, it will be sufficient for our purposes to 
work with g gi ven by Eq. ( 23 ) where g^^^") are defined in 
Eqs. (31 ), (32 ), and the distribution function h presented 



in Eq. (29|. 



IV. CALCULATION OF THE ELECTRIC 
CURRENT 

A. Fluctuation corrections: derivation 

Before studying the fluctuation corrections, we first 
show how to obtain Drude conductivity from the formal- 
ism. Input are the normal-metal solution of the Usadel 
equation: g^ = —g^ = and the distribution function 
in the presence of the electric field, Eq. (29). This gives. 



according to Eq. (22), the electric current: 



(64) 



This results in the Drude formula = 2ve^D. 

Now we turn to the calculation of the fluctuation cor- 



rections. Starting with expression (22), we substitute for 



g the parametrization ( 23 1 and obtain the following con- 
tributions to the current 



.Uc) 



(65) 



Here, all terms besides j^"-' depend on the realization of 
the superconducting order parameter A and have to be 
averaged using the order parameter correlation functions 
(55 ). The fluctuation contributions can be written in the 



following form (hereafter the derivative is with respect to 
the energy argument): 

= 2ne^DEjn'ie)Siy{e){de), (66) 
j^''") = 2ne^DE J H' (e) z9 (e) (de) , (67) 
j(^'=) = 2^eD [ n{e)]^'He)(de). (68) 



The quantities which appear in these expressions are de- 
fined as follows 



die) 



(70) 
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quasiparticle density-of-states: 




and 



FIG. 1: Anomalous Maki- Thompson diagram. 



Jt^ (^) 



(71) 



The rationale behind this decomposition is the follow- 
ing: (i) The function 6i'{e) describes the correction to the 
electronic density of states, see Sv (e) in Eq. ( 72 ) below, 
(ii) The ?5(e)-term has a peculiar analytic structure. In- 
deed, it contains a convolution of /* and /, which upon 
averaging gives rise to a product of two Cooperons of dif- 
ferent analytical structure, and C^, and the imaginary 
part of the fluctuation propagator, ImL^. This allows to 
identify this term with the anomalous Maki- Thompson 
contribution. For an illustration of this point, we refer 
to Fig. [l] (iii) The ji''''-term can be interpreted as the 
fluctuating supercurrent density. This term is the result 
of the expansion in the electric field of the fermionic dis- 
tribution function he entering either the combination V 
(see Eq. (34)) or the order parameter correlation func- 
tion £ (see Eqs. (60) and (63)). The former contribution 



is purely quantum, while the latter comprises both quan- 
tum and classical parts, which are of different importance 
in the different regions of the phase diagram. 

We note that the decomposition (i) - (iii) is very dif- 
ferent from the conventional classification based on the 
diagrams in the Matsubara technique. The difference is 
related to two main points: a) in the traditional tech- 
nique a response to a time-dependent vector potential is 
calculated and b) in the present method there is no need 
for an analytic continuation. 

It is obvious from Eqs. (66) and (67 1 that and 

(an) 
{sc) 



contribute only to the longitudinal current, while 



contributes to the transverse current as well. In 
this context it should be kept in mind that in the Us- 
adel equation, which was used as a starting point for our 
calculation, the Lorentz force acting on the quasiparticle 
was neglected. 

To proceed further, we substitute the expressions for 
/,/*,/ and / in the LL basis (cf. Eq. (40)) into the 
expressions above and average them with respect to or- 
der parameter fluctuations. The quantities 6u(e) and i^(e) 
are equilibrium properties of the system and are indepen- 
dent of the electric field, and that is why their calculation 
is relatively simple. Let us start with the DOS correc- 
tion which can be understood as a renormalization of the 



Si, (e) = Im f {duj) Cl (2e - cj) 

X [Lf(a.)+i»H (6 -..)]. (72) 

Here, v ~ is the number of states per unit area 

of a LL. This factor appears with each summation over 
LLs. In the continuous limit v'Y^^ — ;> and the above 
expression becomes identical to the one in Eq. (372) in 
the review by Kamenev and LevchenkcPH. Note that 
J Sv (e) dt = 0. This is because the interaction cannot 
change the total number of single-particle states, but just 
redistributes them. 

Turning to the anomalous MT correction, we find that 
it is due to a real process. Indeed, i?(e) can be presented 
in the following form: 



ni \ ''"oiit.n(^) 



(73) 



where r^^j „ is the partial in) out-scattering rate for 
quasiparticles arising due to the decay of superconduct- 
ing fluctuations'*^: 



i.„(e) = 2 / (d^) ReC„ (2e - c.) 



xlmL^ (w) {B (uj) uj)] .(74) 

The discussed correction disappears at zero temperature. 
This makes it essentially different from the DOS correc- 
tion which exists down to zero temperature. The sign 
of the anomalous MT correction is always positive. It is 
closely related to weak antilocalization, and can be in- 
terpreted as an interference effect in the singlet Cooper 
channel, enhanced by coherent scattering on the fluctu- 
ating order parameter. 

Next, we turn to the calculation of the supercurrent 
which is more complicated because non-equilibrium 
terms in the fluctuation propagators have to be taken 
into account. The calculation gives: 

jt'W = ^^^E /(da^)(^i+l){A,„+i(^,6)}_(75) 

n 

X {A„,„+i(w,e) - A„,„(a;,e)}^ . (76) 

In these equations, the notation = X ±X \s intro- 

duced, where X is obtained from X by the substitution 
n ■(-> n + 1. The functions Am„(e,cj) are defined in the 
Appendix [B| 

The next step is to substitute 5v (e) , d (e) and ja^ (e) 
into the expressions (66)- (68) and to perform the inte- 
gration in e. The results of these integrations can be 
expressed in terms of £„: 
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= -2e^DvY^ I {du 

n 



B 



.,Imf„Re£'4 



2 



Sa'f'^ = -2e^Dn-^vY^ J (duj) {n + 1) {Bu, + B' v,) , 



(77) 
(78) 
(79) 



where i =||,_L. For the longitudinal (||) conductivity, 

7.11 = Re [if„<i^L«+i] , (80) 
I'll = 2RcKr,lm[£„+Sn+i]lm[L^L^+,] 

+Imif„Im - L^] (81) 



with Kn{uj) = ^^_^_i{uj) — For the transversal (_L) 

conductivity (assuming negatively charged carriers e < 
for the rest of the paper; otherwise, the sign of the Hall 
conductivity should be reversed), we obtain: 



(82) 



= -2Im(^^+^«+i)Rei^„Re [1^1^+,] ~2n, {im^^ImV^f L^L;^}^~Imif„Re {L^^, + L«) +f^,Re {L^Re^f }^ . 

m 



r 



To conclude, we have derived fluctuation conductivity 
due to electron-electron interactions in the Cooper chan- 
nel in the Gaussian approximation. Equations (77l-(79) 



describe the contribution of superconducting fluctuations 
to the conductivity everywhere in the {B, T) phase dia- 
gram (outside the regime of strong fluctuations close to 
the transition) . In the rest of the paper we discuss differ- 
ent limiting cases and elaborate on asymptotics of these 
general formulas. 



B. Discussion: Longitudinal conductivity 

At the end of the previous section, we provided general 
formulas for the fluctuation corrections to conductivity. 
In certain asymptotic regions of the phase diagram they 
are amenable to an analytic treatment. Following this 
route, we are able to compare our results to the previous 
studies. The derived formulas can also be subjected to a 
numerical analysis, which allows to find the corrections 
in the entire normal part of the phase diagram. 

We will discuss the following asymptotic regions in the 
phase diagram: The vicinities of the classical (I) and 
quantum (II) transition points, the region of high tem- 
peratures and small magnetic fields (III) and the region 
of high magnetic fields and low temperatures (IV) . The 
corresponding regions are indicated on the phase diagram 
displayed in Fig. [2j By means of a numerical evaluation. 



we locate the line which describes the transition from 
positive to negative magnetoresistance (9_b(t = 0), and 
the line which characterizes the change of the tempera- 
ture dependence of the total correction drcr = 0. 



GL region (I) 



In this region, '^Cy'*'^'' and Sa^""' are the most impor- 
tant. Since the leading contribution comes from small 
bosonic momenta and frequencies {iu^Dq^ T — Tc), in 
order to extract the result, one should expand the equi- 
librium propagator in cj/T and en/T: 



(an 



8T 



'GL 



where 



tgl 



BTlnT/Tc 



(84) 



(85) 



In this section we assume ^ tql and neglect in 
the fluctuation propagator. Substit utin g the exp ression 



for the propagators Ln^^^ to Eqs. (|78| and (|79|), inte- 



grating in frequency (only the term proportional to B' 
contributes), and performing the summation over the LL 
index, we obtain: 



5CT 



{an) 



—Ttgl 

TT 



( ^ + 5 



(86) 
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FIG. 2: Phase Diagram for the correction to the longitudinal 
conductivity 5(Jxx- The corresponding equations are written 
in the text. 



and 
5a 



(sc) _ 2e2 



with 



-1 - 2stp (s) + 2sil;[ -+s 



{QcTglT 



(87) 



These results are in agreement with existing calculations. 
In particular, Sai^'^^ was obtained phenomenologically by 
Abrahams et al.SS, and the Maki-Thompson contribution 
was discussed for finite magnetic fields in Ref. W, Note, 
that the parameter s divides the region (I) into two parts 
with a distinct behavior. The zero-field limit is recovered 
for s ^ 1: 

^ ^TTGLHr^/TGL), Sa\''^ = ^Ttgl. (89) 

In the absence of a magnetic field, the importance of the 
anomalous MT correction, (5(Tj|""\ in comparison with 

(5(7 Ij''^' is determined by the ratio t^/tgl- Indeed, the 
MT term diverges in the absence of dephasing, — >■ 
oo, and becomes comparable to the AL correction when 
'''4> ^ tgl- As the ratio decreases further, the relative 
importance of the MT correction diminishes. 

For completeness, let us discuss the DOS correction 
in region (I). In the vicinity of the critical temperature, 
g^{dos) -^^gakly (only logarithmically) singular. The 
reason is that interactions preserve the total density of 



states, and the integration with H' in Eq. (661 is (com- 



paratively) wide: e < T k Tc. For zero magnetic field 
one gets: 



5(7 



(dos) 



7C(3)e2 



In Tt, 



GL- 



(90) 



A contribution of the same form originates also from 
the anomalous MT correction as a subleading term, with 
a numerical coefRcient —14 instead of —7. It is instruc- 
tive to perform a comparison with the previously known 
result in this region. For that, one should sum all terms 



of the kind 5cr 



, C(3) 



IhTtgl- In the diagrammatic 



calculation,!^ one obtains the coefRcient c = — 14 as the 
combined contribution of all diagrams with a horizontal 
interaction line. Those diagrams taken together are of- 
ten referred to as the DOS-type corrections. In addition, 
regular MT, AL and anomalous MT diagrams come with 
the coefficients c ~ —7, c = 14 and c = —14, corre- 
spondingly. One can see that only after summation of all 
logarithmic terms of this kind, the results of the two ap- 
proaches coincide, and one obtains in both cases a total 
numerical coefficient Ctot = —21. 



We would like to stress that according to Eq. ( 66 ) it is 



the contribution (5(Ty''°'*' rather than the sum of all hori- 
zontal diagrams that should be associated with the sup- 
pression of the single-particle density of states. 



2. Quantum critical point (11) 

In the vicinity of the transition line, for 

h={B- B^{T))/B^ < 1, (91) 

the most singular contribution conies from the lowest LL, 
n ~ Q. For small temperatures in the vicinity of the 
Quantum Critical Point (QCP), when 



t ^ T/T, < 1, 



(92) 



we can simplify the inverse fluctuation propagator using 
the asymptotic formula for the Digamma function: 



(93) 



£^{u}) = -/i-ln(2n + l) - In 1- 



In this region, the role of in the fluctuation propa- 
gator is mostly to shift the critical magnetic field. We 
will assume that this shift has already been performed. 
Besides, it is natural to neglect in the Cooperon, be- 
cause in the vicinity of the critical point the Cooperon is 
not singular and I/t^ has to compete with fi^. Substi- 



tuting the expression for i£'„(a;) into Eqs. (77l-(79) and 
expanding the propagators in to file, the contributions of 
all three terms can be written in the form 



5a^ 



(0 



o2 r 



.W/„(t,/i)-f /3«/^(i,/i) 



(94) 
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with the numerical coefficients 

^{dos) ^ _j_ 



=0, = ^, (95) 
3 

(dos) ^ p{an) ^ 2, jSi'") = ^. (96) 

3 



Here 



He 



ujB(w)duj 



Evaluating these integrals, we obtain: 
(t, h) 
Ip {t, h) 

with 



w2 + (hVlj2Y 



r 1 

n. 2r 

rT/;' (r) 1 



1 h 



(97) 
(98) 

(99) 



Note that when all the contributions are summed up, we 
get a = — |, /3 = |, and our result reproduces the one 
obtained by Galitski and LarkirP. 

The region of the phase diagram in the vicinity of the 
QCP can further be subdivided into classical and quan- 
tum regions, depending on the ratio of the parameters h 
and t. The superconducting fluctuations contribute ei- 
ther as classically populated modes or through virtual 
transitions. In the quantum region t <^ h the occupation 
number of the lowest LL of the collective mode is small, 
and we obtain 



Sail 



2e2 1 



{t < h). 



(100) 



In the classical region t ^ h, the occupation number 
is large and the correction changes its character. As a 
result, it becomes positive: 



6(7n 



26^7 t 



(t > /i). 



(101) 



3. High temperatures (III) and high magnetic fields (IV) 

In these regions the dominant contributions come from 
high LLs and, hence, the summation in the LL index can 
be replaced by an integration. At the same time, the full 
dependence of the fluctuation propagators on the bosonic 
frequency should be kept, because the leading contribu- 
tion comes from a long double logarithmic integration. 

Let us first discuss the region (III). We will perform 
the calculation in the limit of ln(r/rc) ^ 1- We start 



with the analysis of Sa 



(dos) 



It has a very slow temper- 



ature dependence due to the long integration in energy, 
which has to be cut off at oj, e ~ t~^, where the diffusive 
approximation breaks down. In view of this fact, only the 



term proportional to B (rather then B') gives the leading 
contribution, and we can write 



Sa 



(dos) 



e 

4^ 



B (w) Im [i^ (uj) ?A^"(w)] dujde 

B{uj)dl^{h + '-^]dude 



47r2 



Im 



lnT/T, + V 



e — iuj ^ 
47rT ) 



(102) 



This integral is logarithmically divergent. As a result, we 
obtain: 



" 27r2 InT/Tc 



(103) 



This correction is similar to the Altshuler-Aronov cor- 
rections, but with a scale-dependent coupling constant. 
This result was first derived by Altshuler et At 
very large temperatures (InT/Tc 3> 1) this term dom- 
inates the total correction. In the case of a repulsive 
interaction, it become^^ In In J- . 

Let us turn to The term proportional to B' 

is again small, O (In^^ (T/Tc)). Another term, which is 
proportional to B, is more important: 



izdz 



d?/cothfV' (e)V'" (e) 



256^5 [lnT/Te + V(£), 

(104) 

where e = ^ + ^j^- We first calculate the y integral 
neglecting y in the denominator. Since only y ^ 1 con- 
tribute to the leading term, we can substitute coth | ^■ 
signy. This leads to 



647r4 



zdz [i^'i'^ + i,)]' 
[lnr/T, + ,^(i + ^)] 



2 • 



(105) 



The remaining integral comes from 1 ^ z and can be 
calculated to give: 



5(7 



(sc) 



1 



47r2 InT/Tc 



(106) 



We note, however, that the same term originates from the 



subleading contribution to 5(7}^'^'"'\ but with a different 



numerical coefficient 



In 2-1 
27r2 



Thus, different contributions 



of the kind 0{\n~^ T/T^) do not cancel each other. 

Let us now turn to (5cr|™\ In the continuous limit, 
vJ2n ~^ J2qy Eq. (78) reproduces the known resuliP^. 



In the limit of hiT/Tc ^ 1, it can be further simplified 
to: 



Sa 



(an) 



1 



with 



M{z) 



M{z)dz 

167r2 In'YJr, Jo z + 1/ (Tr^) 
- dy[^{l 



(107) 



z-ry\ 
4tt ) 



4, 



sinh^ {y/2) 



(108) 
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Although this term is formally ©(In"^ T/T^), it can still 
be essential due to the logarithmic divergence at small 
momenta as it can be seen from Eq. (107). With loga- 



rithmic accuracy, we can calculate it as follows 



5a 



(an) 



1 



M{{))dz 



167r2 T/T, Jo z + l/ (Tt^) ' 
As a result, we get: 



(5(7; 



(an) 



In Tr^ 
12 In^ T/Tc 



(109) 



(110) 



One should keep in mind, however, that itself depends 
on T. In this region, the anomalous Maki-Thompson 
correction was considered by several authors, who all ob- 
tained the same functional form but with different nu- 
merical coefficients^"^ We believe this discrepancy 
is due to different approximations used for the calcula- 
tion of Af(0). 

For high magnetic fields (region (IV)), the situation is 
to some extent analogous to region (III) with the main 
difference that the anomalous MT term does not con- 
tribute as it is suppressed at small temperature. The 
dominant corrections originate from Sa'^^'^'' and 5a^'^°^\ 
and the leading contributions are those which are pro- 
portional to S « signcj. To proceed, we write the equi- 
librium propagator in its zero-temperature form: 



L«(^) = In- 



(T^ 0). 



(Ill) 



After the frequency integration, we find that Sa^^°^' takes 
the following form: 



(112) 



with the logarithmic integral function li(2;) = dt/lnt. 
This sum is logarithmically divergent at the upper limit 
and has to be cut off when the diffusion approximation 
breaks down, that is at n ^ Nmax 3> 1 with N^ax — 
j-jr^- Under these conditions, the sum is dominated by 
large n and can be found to equal 



5(7 



{dos 



In 1/tT, 
27r2 "^InB/B^ 



(113) 



This concludes our discussion of the regions (I-IV) in 
the phase diagram; the corresponding asymptotic expres- 
sions are referenced in Fig. [2j 

The results we obtained differ from those given in 
Ref. |23l This follows from a comparison of the asymp- 
totic behavior in several regions. The most drastic differ- 
ence, however, concerns the temperature dependence of 
the resistance for magnetic fields B > Be- The authors 
of Ref. [23| claimed that for small temperatures T ^ Tc 
the resistance first increases with increasing T and starts 
to diminish at T/Tc > {B - Bc)/Bc. As follows from 



0^ 




0.25 



0.5 

T/T, 



0.75 



FIG. 3: Resistance as a function of temperature for magnetic 
fields B/ Be = 0.9, 1.05, 1.1, 1.3. The sample parameters are 
Rd = hkO. and T^r = 10"^ 



our asymptotic expressions presented in Eqs. ( |94[ ) and 
from the result of the numerical calculation shown in 
Figs. [2]and|3j the situation is opposite. At a fixed mag- 
netic field, the resistance decreases as the temperature 
increases from zero until the line 9t(t = is crossed. 
Then the resistance starts to grow. 



C. Discussion: Hall conductivity 

We proceed with the discussion of the results for the 
transverse conductivity presented in Eq. ( 79 ) . These 



expressions represent only those contributions to 5a ±, 
which describe a deflection of the supercurrent. In prin- 
ciple other contributions exist, in which quasiparticles 
are deflected in the transverse direction by the Lorentz 
force. These contributions are not included in the ap- 
proximation we apply here. The terms not accounted for 
by Eq. ( 79 ) include the contribution due to the anoma- 
lous MT process, discussed by Fukuyama et all^and the 
contribution 5a'^^°^\ recently discovered diagrammati- 
cally by Michaeli et al.l^, which is reminiscent of the 
density of state suppression. They are related to the cor- 
responding corrections to the longitudinal conductivity 
as follows: 



5a 
5a 



(an) 



(dos) 



Note, that 5a^^^'' and Sa^^"""' cancel each other in the 
expression for the Hall resistivity p^y = —<^xy/{<^xx + 
'^xy) ~ ~^xy/<ytx- III contrast, the DOS-corrections give 
a finite contribution to pxy 

Let us discuss the contribution to Hall conductivity 
that arises due to the deflection of the fluctuating su- 
percurrent. In order to calculate it, it is enough to mod- 



r, jr (an) 
-'2uJcT5a]^ , 

r (dos) 

5a „ '. 

2 II 

(an) 



(114) 
(115) 
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if y th e superconducting fluctuation propagator according 



t( 



J28l 



r-1 



r-i 



(116) 



As a consequence of the additional term, the supercon- 
ducting propagators lose their particle/hole symmetry, 
i.e., the relation L^{—uj) ~ L'^^lu) no longer holds. In 
the framework of the BCS theory, the asymmetry pa- 
rameter q can be related to the energy dependence of 
the density of states at the Fermi level: <; — —^^j^f- 
or, equivalentljEI!^ to the variation of with the chem- 
ical potential: <; — —j^^t^- In the simple model of 
3D electrons with a quadratic spectrum, one has i'(e) « 
1^0(1 + e/2eF) and <: = -l/(4eFA). For A < 1 the con- 
tributions arising from (5(7^'^'' are parametrically larger 

than those arising from (5(7^°'''' and Sa^^"\ In our calcu- 
lation of the Hall conductivity, we work in the framework 
of the quasiclassical approach, using, however, Eq. (116) 



for the propagators L^^a) ■ This is a consistent procedure 
that allows to obtain all contributions to the transverse 
current proportional to the large parameter 1/A. 

In region (I) after expansion in Vlc{n + l/2)/47rT and 
lo/A-kT the correction ^cr^'^'' takes the form 



(117) 



where 



V' Q + - V' (1 + s) - s%Ij' (1 + s) 



In this region, the Hall effect can be considered phe- 
nomenologically: the same expression ( 117[ ) was obtained 
by Aronov and Rapoport^ (with a different coefficient, 
it has later been corrected by Aronov et al.!^) on the ba- 
sis of the time dependent Ginzburg-Landau theory. For 
s 3> 1, when quantization of the LLs for the supercon- 
ducting fluctuations is negligible, the expression (117) 
becomes!^ 



5a 



(sc) 



96 



T 



(118) 



The region of applicability of the Eq. (1171 is in fact 



very narrow, and already for T > l.OlTc one should not 
expand the full expression for (Jcr^^' in flc(n+ l/2)/47rT 
to get an accurate result. The corresponding formula has 
been given in Ref. [231 



da'l 



(sc) 



(0) 



In region (II), we can limit ourselves to the lowest LL 
and follow the same route as in the calculation of the 
longitudinal conductivity. This gives for the quantum 
regime: 



Sa 



(sc) 



37r2 



In 



1 



(120) 




FIG. 4: Resistance as a function of magnetic field for temper- 
atures T/Tc = 0.03, 0.1, 0.35. Inset: the zoomed region of 
the approximate crossing for T/Tc = 0.15 — 0.3. The sample 
parameters are Rd = 5kQ, and TcT — 10^^. 



and for the classical regime: 



5a 



(sc) 



2e2 



(121) 



Note, that in this region 5af^°'^'^ and Jcr^"-' exhibit the 

(sc) 

same singular behavior as <5crj_ . We do not provide 
the corresponding expressions, since they follow straight- 
forwardly from Eqs. (114) and (115), together with 
Eq. (Igil). 



A more detailed discussion of the corrections to the 
Hall conductivity due to superconducting fluctuations is 
presented in a separate publication, Ref. 



V. CONCLUSION 

We considered homogeneously disordered films above 
the superconducting transition T > Tc{B) and calculated 
corrections to longitudinal as well as transversal conduc- 
tivities. Our results are presented by equations ( 77 )- ( 79 1 . 



We analyzed the asymptotic behavior of these corrections 
in different regions of the phase diagram and provided a 
comparison with previously published results. 

Our results for the Hall effect have recently been used 
in the description of experimental data by Breznay et. 
aP3. The results for the longitudinal conductivity, Eqs. 



119) (77)-(79), can also be useful for the analysis of experi- 



ments. They allow for a complete numerical evaluation 
of the fluctuation corrections to conductivity without 
any additional approximation, e.g., the lowest Landau 
level approximation. Exemplary results are presented in 
Figs. kJI UJfor the resistivity R = {R~j^ + 5a)~^ as a func- 
tion of magnetic field and temperature. A similar behav- 
ior of the resistance was observed in the experiment of 
Baturina et In Ref. 151 the authors presented a fit 
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to the measured data that was based on the asymptotic 
expressions ( 94 ) derived in Ref . 14J and reproduced in our 
work based on a different method. We note, however, 
that although these expressions provide a good approxi- 
mation in the vicinity of the QCP, their region of vahdity 
does not extend up to the relatively large temperatures 
and magnetic fields that were considered in the exper- 
iment (up to O.SSTc and up to 5Bc, correspondingly). 
When fitting this data, the more precise Eqs. ([77|-(79) 
should, therefore, be used. 

According to the results presented in this work, the re- 
sistance curves drawn as a function of the magnetic field 
exhibit an approximate crossing point for a finite interval 
of temperatures, as demonstrated in Fig. [3) As can be 
seen from this picture, the curves do not literally cross 
in a single point, but deviations from this ideal behav- 
ior are small. The existence of this approximate crossing 
point is a consequence of a relatively wide minimum in 
the R{T) curve for B = 1.053^ as shown in Fig.[3| This 
type of behavior has been observed in several systems; 
see e.g. Fig. 4 in Ref. [52l However, in these experiments 
the curves continue to cross even at the smallest temper- 
atures, while we did not find this kind of behavior from 
the Gaussian corrections to conductivity. This could be 
related to the fact that for such low temperatures the 
proximity to the QCP becomes of crucial importance, 
and the present theory is not sufficient because 1) it does 
not account for the effect of non-Gaussian fiuctuations 
and 2) does not take into account the smearing of the 
transition by disordeJ^^'^, which is usually observed in 
this region (see Fig. 2 in Ref. |S5]as an example). 

To conclude, we have developed an approach to the cal- 
culation of fluctuation conductivity based on the Usadel 
equation and valid for both the classical as well as the 
quantum fluctuation regime for arbitrary magnetic flelds. 
This approach is more physically transparent than con- 
ventional perturbation theory based on the Kubo formula 
and provides a bridge between the phenomenological the- 
ory and microscopies. We believe that it may flnd appli- 
cations in studies of fluctuation effects out of equilibrium 
or in hybrid superconductor/normal metal structures. 
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Appendix A: Collision integrals 

In this appendix we discuss the quantum components 
jZ,w -j-j^g Qi-een's function g^. We parameterize them 



9 



zi 

Z2 



9 



w 



wi 

W2 



and get the following equations: 
with 



(Al) 



(A2) 



DV^ - dt, - 5t,, = DV^ + dt, + dt,. (A3) 



The collision integrals 2 are given by 

/f = z(A,.r-/.A;) 

/f = z (a: . / - f . A,) 



and collision integrals /,. 
ration is motivated below) 



w _ tW 



i,coll 



tW 

i^neq 



(A4) 

by (this sepa- 



tW 

tW 

'^2,coll 



«(/• Ji- Ji-f), 
= *(/* • J2 - J2 • /), 



w 

neq 



1. 



w 



f'-jh-r + f'-r-je, 

.2,ne, - ^Jh ■ Z2 ■ Jh + Jh ■ r ■ f + r ■ Je ' /' + 

r-je-f+r-f-jh- (A5) 

For convenience, we deflned {je.h — ±V/ie./i): 
Ji = Al-Al-K + hh-Al-hh-A*-h,, 



Ji 
J2 
J2 



Ac ■ hh + he ■ Ac 

Ac ■ hh + he ■ Ac - IOC ■ 

Al-hc + hh ■ Al - hh ■ A* ■ he. (A6) 



While (/^) — due to causafitjJSS^ the collision inte- 
gral does not vanish identically after averaging. Nev- 
ertheless, its expansion in the electric fleld can be shown 
to start from E^. First, we note that ij^eq should be 
related to the production of the heat. Indeed, (jYneql 
is proportional to the Drude result for the electric cur- 
rent ie,h- Next, observe that the terms in (l^neql which 
are only linear in 2e.,h are further multiplied by averages 
which include the spatial gradients of / and vanish in the 
absence of an electric field, when the system is isotropic. 
Hence, (/^e,) = ©(E^). There is stiU another term, 
lYcoii- ^'^^ E = it corresponds to the collision integral 
due to Cooper interactions, which enters the kinetic equa- 
tion and was calculated by Reizei"*^. Let us just note, 
that if the only source of non-homogeneity is a spatially 
varying electric potential (as it is in our case), then the 
collision integral, written in terms of the gauge invariant 
particle/hole energies should be independent of the spa- 
tial coordinates. As such, it cannot depend on the electric 
field itself, which is a vector, but only on E^. This is sum- 
marized by the equation: ^^^j — IcoIi(Ej^, e T e0(x)). 
Since for E = it vanishes (provided the electronic dis- 
tribution function T-l is thermal) and depends only on E^ , 
it should be disregarded for the calculations in the linear 
response. 
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Appendix B: Calculation of the supercurrent 

Here we present more details of the calculation of 
ja \e)- We start with expression ( 71 ). After substituting 



the solution for / and averaging in A we get: 

ji'^ (e) = \^EY. j idu^)Ic.,rr.nA^mni^, e) (Bl) 
mn 

Here Ia,mn represents the result of integration in the mo- 
mentum quantum number: 

la.rnn = 2j J {dp)lm (^Vmp(r) Vq ^^^(r)^ {np\ X \mp) 

(B2) 

and Ajnn — X]fe Amn has several contributions, which 
arise from different ways to expand propagators or 
bosonic/fermionic distribution functions in the electric 
field. The next step is to calculate integral (B2): taking 



into account {n,p\ x \m,p) = Xnm + pl'sSnm, where a;„, 
are matrix elements, calculated with Xnix), we obtain: 



Iy{m, n) 



2i 

'■B 



We also take into account: 



Xmn = + I6m.,n+1 + VnSm,n-l) 



V2' 



V2I1 



(B3) 



(B4) 



(B5) 



m^n— 1 ) (B6) 



and obtain the result, presented in (75), (76) 
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